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Abstract 

In the framework of the multiconfiguration Dirac-Hartree-Fock method, we investigate the tran- 
sition properties of four excited states in the 2p 5 3s configuration of neutral neon. The electron 
correlation effects are taken into account systematically by using the active space approach. The 
effect of higher-order correlation on fine structures is shown. We also study the influence of the 
Breit interaction and find that it reduces the oscillator strength of the 3 P^ — 1 Sq transition by 
17%. It turns out that the inclusion of the Breit interaction is essential even for such a light 
atomic system. Our ab initio calculated line strengths, oscillator strengths and transition rates are 
compared with other theoretical values and experimental measurements. Good agreement is found 
except for the 3 P$ — 1 Sq M2 transition for which discrepancies of around 15% between theories 
and experiments remain. In addition, the impact of hyperfine interactions on the lifetimes of the 
3 Pq and 3 P£ metastable states is investigated for the 21 Ne isotope (1=3/2). We find that hyperfine 
interactions reduce the lifetimes drastically. For the 3 Pq state the lifetime is decreased by a factor 
of 630. 



PACS numbers: 31.15.A-, 31.15.ag, 31.30.J-, 32.10.Fn 



I. INTRODUCTION 



Lifetimes of states in the first excited configuration 2p 5 3s for neutral neon are important, 
not only because of their potential applications [l-7] relevant to plasma diagnosis, laser 
techniques, and the interpretation of astrophysical data, but also for investigating electron 
correlation and relativistic effects or testing many-body theories of atomic physics |8- 14] . 



A great deal of calculations and measurements have been reported for electric dipole (El) 
transition probabilities or corresponding oscillator strengths (gf) between the states of the 
ground 2p 6 and first excited 2p 5 3s configurations of neutral neon. However, a satisfactory 
precision has not been achieved yet. For instance, the gf values of the lower J = 1 level, 
i.e. 3 P° (the LS coupling label is used throughout this paper for convenience), obtained 
by the nonrelativistic wave functions with relativistic corrections in the Breit-Pauli (BP) 
approximation range from 0.0102 to 0.0123 , while the relativistic results are 

larger than 0.0130 jll[ [3]. Unfortunately, the inconsistency cannot be resolved by the 
experimental measurements because of large error bars. The detailed comparisons have 

n n n n 

been recently reviewed by Chan [15], Avgoustoglou [111 ]. Savukov [12] and Zatsarinny [141 ] . 

Another appealing subject is the lifetimes of the two metastable S P£ and 3 Pq levels in 
the 2p 5 3s configuration. For isotopes without nuclear spin /, the magnetic quadrupole (M2) 
transition to the ground state is the dominant single-photon decay channel for the state, 
while the 3 Pq level can decay through the magnetic dipole (Ml) or electric quadrupole (E2) 
transition to 3 P° 2 lower states. In 1972, Van Dyck, Johnson, and Shugart measured the 
composite lifetime of the metastable rare-gas atoms in these two states using the time- 
of-flight technique 16j. The experiment sets a lower limit for the lifetime, and the value 
is 0.8 s for Ne. Recently, Zinner determined the lifetime of the 3 P 2 ° state by measuring 
the decay in fluorescence of an ensemble of 20 Ne atoms trapped in a magneto-optical trap 
(MOT) 17]. It is worth noting, however, that the latest experimental result r = 14.73(14) s 
considerably differs from the earliest theoretical results^ r = 24.4 s by Small- Warren and 
Chow Chiu [18] and r = 29 s by Fielder, 
recent calculations, that is, 19.8 s by Beck 



r. et al. [8|, [9J. Also it does not agree with 



19] with relativistic configuration interaction 



method, 18.9 s obtained by Desclaux et al. (cited in [17] ) and Dong et al. [13] using the 
multiconfiguration Dirac-Hartree-Fock (MCDHF) method and 16.9 s by Froese Fischer and 
Tachiev [5] with the multiconfiguration Hartree-Fock (MCHF) method including relativistic 



corrections in the BP approximation. 

On the other hand, for isotopes having a non-zero nuclear spin, issues become complicated 
since a new decay channel is opened by hyperfine interactions. This transition is referred to 
as a hyperfine induced transition (HIT) or hyperfine quenching decay mode [20j]. Owing to 



their peculiarity, HITs have attracted much attention during the last several years [2lN25|. 
stimulating us to further predict the rates for the 3 Pq an d 3 P$ metastable states of the 21 Ne 
isotope. 

re transition properties of states in 



26| based on the MCDHF method 



In this work, we perform large-scale calculations of t 
the 2p 5 3s configuration using the GRASP2K package z 
which allows one to take electron correlation and relativity into account on the same foot. 
The active space approach is adopted to monitor the convergence of the physical quantities 
concerned. The importance of the Breit interaction for an accurate determination of the 
lifetimes is studied. We report the lifetime of these states for abundant isotopes with respect 
to important decay channels including hyperfine induced transitions. 

II. THEORETICAL METHOD AND COMPUTATIONAL MODEL 
A. MCDHF method 

The multiconfiguration Dirac-Hartree-Fock method is written up in the monograph by 



Grant 22J and we here just give a brief description of the method. Starting from the Dirac- 



Coulomb Hamiltonian 

h dc + & ~ l )° 2 + O + XJ Vry, (!) 

i i>j 

where V N is the monopole part of the electron-nucleus Coulomb interaction, the atomic state 
functions (ASFs) describing different fine structure levels are obtained as linear combinations 
of symmetry adapted configuration state functions (CSFs) with same parity P, angular 
momentum J and its Mj component along ^-direction 

JV 

V(PJMj) = Y^cM^PJMj). (2) 

i=i 

In Eq. (j2j) Cj is the mixing coefficient and jj denotes other appropriate labeling of the 
configuration state function, for example orbital occupation numbers and coupling tree. 



The configuration state functions are built from products of one-electron Dirac orbitals. In 
the self-consistent field (SCF) procedure both the radial parts of the Dirac orbitals and the 
expansion coefficients are optimized to minimize the energies concerned. Calculations can 
be performed for a single level, but also for a portion of a spectrum in an extended optimal 
level (EOL) scheme where optimization is applied on a weighted sum of energies. The Breit 
interaction 
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can be further included in subsequent relativistic configuration interaction (RCI) computa- 
tions. 

Once the atomic state functions have been obtained atomic parameters are evaluated in 
terms of reduced matrix elements of the corresponding tensor operator 

( ^(PJ) ||0 (A) || ty(P'J') ). (4) 

For the transition, the tensor operator is a multipole radiation field operator. The 
superscript designates the type of multipole: A = 1 for electric multipoles and A = for 
magnetic multipoles. This expectation value reduces to a sum over reduced matrix elements 
between CSFs by substituting the ASF expansions (J2]). Using Racah algebra these reduced 
matrix elements, in turn, are expressed as a weighted sum over radial integrals involving the 
radial relativistic one-electron orbitals. 

The restriction from Racah algebra that ASFs are built from the same orthogonal radial 
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29|. As a result, 



orbital set can be relaxed by the biorthogonal transformation technique 
reduced matrix elements between two atomic state functions described by independently 
optimized orbital sets can be calculated using standard techniques. 



B. Computational model 

In the framework of the MCDHF method, the building of the configuration space is 
pivotal not only for capturing the electron correlation effect efficiently, but also for circum- 
venting the convergence problem that one frequently encounters in SCF calculations. In this 
work, we use the active space (AS) approach to generate the configuration list from the ref- 
erence configuration set. The reference set is initially made up of (near-) degenerate reference 
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configurations and can be augmented by important CSFs for considering the higher-order 



correlation effects 3CH32I]. We name in this paper the initial set as MR^ and the latter 
MRW. More generally, the reference set MR can be divided into several subsets for ex- 
plaining correlation effects between specific electron pairs. According to the perturbation 
theory, the first-order correction of ASFs is expressed as a linear combination of CSFs that 
are obtained by replacing one or two occupied orbitals of the reference configurations in 
MR' ' with active orbitals 3s|. The set of active orbitals is enlarged systematically, which 
makes it possible to monitor the convergence of the physical quantities under investigation. 
Higher-order correlation corrections are more difficult to deal with since the number of CSFs 
grows rapidly and easily goes beyond the capability of even a large computer system. Yet 
most of CSFs actually make fractional contributions to ASFs. The key point in this step is 
to define the MR' 1 ' appropriately. In general, significant CSFs in first-order correction are 
added to the MR' ' to form the MR' 1 ' set. The configuration space is further expanded by 
single(S) and double(D) replacements for orbitals of CSFs belonging to MR 1 - 1 ) with the ones 
appearing in a given active set. 

For the case of neon, Lee et al. pointed out that higher-order correlations of the L 
shell are significant for the ground state [341 ] . Afterwards, Dong et al. also showed by 
MCDHF calculations that CSFs generated from the 2s 2 2p 4 3p 2 configuration improves 



the accuracy of the transition rates to a great extent (l3j. As a result, we choose 
{{2s 2 2p 6 ;2s 2 2p 4 3p 2 ;2s 2 2p 5 3p} ; {ls 2 2s 2 2p 6 }} as the MR(°' set for the ground state and 
{{2s 2 2p 5 3s} ; {ls 2 2s 2 2p 5 3s}} for the four lowest excited states, respectively. The first sub- 
sets in MR' ' aim at accounting for the outer electron correlations and the second for 
correlations involving the Is core. The construction of the configuration space is presented 
in table [H As can be seen from this table, these correlation models are marked with n/SD 
where n and I, if appearing, designate, respectively, the maximum principal and orbital an- 
gular momentum quantum numbers of the active orbitals. The core correlation involving Is 
electrons (labeled as "CC") is taken into account by allowing SD excitations from the Is core 
to the largest active set. To incorporate the residual higher-order correlations of outer shells, 
the {2s2p b 3s3p ; 2s 2 2p A 3d 2 ; 2s2p 6 3s} configurations are added to the first subset of MR' ' 
for the ground state and {2s 2 2p 3 3s3p 2 ; 2s2p 5 3s3d ; 2s 2 2p 3 3s3d 2 ; 2s 2 2p 4 3s3p ; 2s 2 2p 5 3d} for 
the excited states to set up MR' 1 ' . The SD excitations up to n = 4 are based on the MR' 1 ' 
and the CSFs are appended to the CC model to form the final configuration spaces (marked 



with MR). It is worth noting that the addition of 2s 2 2p 5 3p ; 2s2p 6 3s ; 2s 2 2p 4 3s3p ; 2s 2 2p 5 3d 
configurations in the reference sets is ascribed to the requirement of closing the CSF space 
under de-excitation by the biorthogonal transformation technique 26] . 



In practice, we further eliminate the CSFs that do not interact with reference config- 



urations 
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331 ] in order to raise the calculation efficiency. As can be seen from table [I] 



the number of CSFs of the reduced configuration space is considerably smaller than the 
corresponding full one. These removed CSFs contribute to the atomic properties under 
investigation at higher-order and the quantitative influence can be seen in table HD Using 
9fSD and CC models, we calculate excitation energies and 3 P° — 1 S'o line strength. It is 
found that the impact of removed CSFs on excitation energies between different configu- 
rations or terms is fractional whereas remarkably large for the fine structure splitting. For 
example, the influence reaches around 20% for the 3 P° — 3 P£ fine structure. Comparing 
the El line strengths obtained with the two configuration spaces, we see that the loss of 
CSFs contributes to the weak line by 3% but merely 0.2% for the strong line. 



C. Breit Interaction 



Ynnerman et al. have demonstrated that the Breit interaction plays a key role in the 



spin-forbidden 3 P° — l So transition of low- Z Be-like ions 36]. Avgoustoglou et al. have 
also illustrated the effect of the Breit interaction on the Ne I transition energies |lo|. In 
this subsection, we investigate the Breit interaction effect on transition energies and on the 
3 P° — 1 S line strength. It should be emphasized that the full configuration space must 
be used because the strategy adopted for reducing the number of CSFs does not apply to 
the Dirac-Breit Hamiltonian. As examples, we present results with and without the Breit 
interaction, which are obtained using the full DF, 9fSD and CC configuration models. As 
can be seen from this table, the Breit interaction substantially affects the physical quantities 
concerned. For instance, the impact of the Breit interaction on the line strength for the 
3 P° — 1 S transition reaches about 17%. 



TABLE I. The number of configuration wave functions (CSFs) in various correlation models. J p are 
the total angular momentum (J) and parity (P) of an atomic state. MR stands for the reference 
configuration set, AO for the set of active orbitals. The number of CSFs without reduction is 
presented in the parentheses following the number of the reduced configuration space. * indicates 
that all active orbitals are included. 



NCSF 



MR 


AO 


Model 




J p =0 e 




{2s 2 2p 6 ; 2s 2 2p 4 3p 2 ; 2s 2 2p 5 3p} 




DF 




12(12) 






{3s,3p,3d} 


3SD 




549(728) 






{3*,4s,4p,4d,4f} 


4SD 




3731(6021) 






{3*,4*,5s,5p,5d,5f,5g} 


5SD 




10884(19355) 






{3*,4*,5*,6s,6p,6d,6f,6g,6h} 


6SD 




23166(43967) 






{3* 4*,5*,6*,7s,7p,7f,7f} 


7fSD 




35746(67433) 






{3*,4*,5*,6*,7*,8s,8p,8d,8fj 


8fSD 




51122(96017) 






{3*,4*,5*,6*,7*,8*,9s,9p,9d,9f} 


9fSD 




69294(129719) 




1 1/1 <! 2 2<s 2 2r> 6 l. 


{3s,. . . ,9f} 


CC 




71406(132005) 




\J{2s2p 5 3s3p ; 2s 2 2p 4 3d 2 ; 


{3s,3p,3d,4s,4p,4d,4f} 


MR 




81327(143037) 




2s 2 2p 5 3p;2s2p 6 3s} 


















J p =0° 


J p = 1° 


J p = 2° 


{2s 2 2p 5 3s} 




DF 


1(1) 


2(2) 


1(1) 




{3s,3p,3d} 


3SD 


86(145) 


326(369) 


287(431) 




{3*,4s,4p,4d,4f} 


4SD 


444(866) 


1942(2279) 


1821(2887) 




{3*,4*,5s,5p,5d,5f,5g} 


5SD 


1192(2495) 


5500(6734) 


5327(9027) 




{3*,4*,5*,6s,6p,6d,6f,6g,6h} 


6SD 


2442(5325) 


11600(14639) 


11469(20435) 




{3* 4*,5*,6*,7s,7p,7f,7f} 


7fSD 


3727(8162) 


17846(22332) 


17611(30868) 




{3*,4*,5*,6*,7*,8s,8p,8d,8f} 


8fSD 


5289(11619) 


25468(31683) 


25100(43485) 




{3*,4*,5*,6*,7*,8*,9s,9p,9d,9f} 


9fSD 


7128(15696) 


34466(42692) 


33936(58286) 


U{ls 2 2s 2 2p 5 3s} 


{3s,. . . ,9f} 


CC 


11744(30740) 


59320(83520) 


55901(113950) 


(J{2s 2 2p 3 3s3p 2 ; 2s2p 5 3s3d ; 


{3s,3p,3d,4s,4p,4d,4f} 


MR 


45368(63831) 135830(173967) 


187309(238761) 



2s 2 2p 3 3s3d 2 ; 2s 2 2p 4 3s3p ; 2s 2 2p 5 3d} 



III. RESULTS AND DISCUSSION 



A. Excitation Energies 

As functions of the computational models described in Section III Bl the excitation ener- 
gies are presented in table |TV] The reduced configuration space is used at each step except 
for the last one where the Breit interaction is accounted for with the full configuration list. 
It is found that the correlation between outer electrons is saturated in the 9fSD model. 
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TABLE II. Comparison of transition energies (in cm -1 ) and El line strengths (in a.u.) obtained 
with reduced (r) and full (f) 9fSD and CC configuration spaces. B: Babushkin gauge; C: Coulomb 
gauge. NIST data [3] are presented for reference. 







Excitation energy 




Line strength ( 3 Pf - l So) 


Line strength ^Pf - 1 S ) 


Model 


'So - 3 P? 


3 po 1 po 


3 po 3 po 
r X M) 


3 po 3 po 
r 1 M) 


B 


C 


B 


C 


9fSD(r) 


134838 


1391.95 


475.16 


874.87 


0.03417 


0.03419 


0.3562 


0.3556 


9fSD(f) 


134837 


1397.68 


400.33 


835.33 


0.03512 


0.03508 


0.3555 


0.3550 


CC (r) 


135398 


1381.44 


479.10 


880.12 


0.03488 


0.03410 


0.3524 


0.3440 


CC (f) 


135395 


1387.28 


404.84 


840.76 


0.03583 


0.03499 


0.3517 


0.3435 


NIST 


134459 


1429.43 


359.35 


776.80 











TABLE III. Breit interaction effect on the transition energies and the 3 P° — 1 Sq line strength. 
The full configuration space is used in these calculations. B: Babushkin gauge; C: Coulomb gauge. 
NIST data [35J are presented for references. 





Fine structures (in 


cm 1 ) 


Line strength ( 3 P° — 1 Sq) (in a.u.) 


Model 


3 Pf - x Pf 


3 po 3 po 
r l " M) 


3 po 3 po 
r 2 ~ M) 


B 


C 


DF 


1400.06 


389.16 


820.83 


0.02584 


0.03172 


DF + Breit 


1485.00 


347.73 


752.63 


0.02184 


0.02684 


9fSD 


1397.68 


400.33 


835.33 


0.03512 


0.03508 


9fSD + Breit 


1381.77 


358.59 


767.36 


0.02989 


0.02986 


CC 


1387.28 


404.84 


840.76 


0.03583 


0.03499 


CC + Breit 


1370.59 


362.71 


772.65 


0.03054 


0.02979 


NIST 


1429.43 


359.35 


776.80 







The core correlation and high-order effect make relatively small contributions but signif- 
icant to bring the excitation energies to a satisfactory agreement with the experimental 
values 35]. Compared with other theories, the present excitation energies between ground 
and excited states are better than those obtained by Avgoustoglou et al. [10|, 111] and by 



Savukov et al. 12] with many-body perturbation theory (MBPT), but are not as excellent 
as the MCDHF data of Dong et al. [lsj and MCHF values of Froese Fischer and Tachiev jsj. 



S 



TABLE IV. Excitation energies (in cm 1 ) of 2p 5 3s levels for neutral neon. 



Model 


3 DO 


3po 


3 no 
M) 


lpo 


3 Pf - 


3 no 3 no 


3 no 
r 2 — 


DF 


140733 


141165 


141554 


142565 


1400 


389 


821 


3SD 


138180 


138599 


138986 


139977 


1378 


387 


806 


4SD 


133797 


134203 


134652 


135600 


1397 


449 


854 


5SD 


134420 


134821 


135288 


136207 


1386 


468 


869 


6SD 


134572 


134972 


135446 


136352 


1380 


474 


875 


7fSD 


134494 


134894 


135370 


136277 


1383 


475 


876 


8fSD 


134462 


134861 


135336 


136253 


1392 


475 


874 


ytoJj 


1 Q/1 A QQ 


134838 


135313 


lOOZOU 


1 QOO 


A IK 

4(0 


Q7CT 




1 Q/1007 


135398 


135877 


1 Qfivcn 
loD/oU 


1 QQO 


4 /y 


ooU 


IVlxt 


Io4o4 ( 


134783 


135191 


1 Qfil VQ 


i Qon 


4Uo 


840 


Breit 


1 Q/l Q^fi 


134765 


135127 


1 Qfil A 1 


1 Q7^ 


QR9 


771 
111 






Others 










Avgoustoglou et al. [10] 


134011 


134406 134757 


135570 


1164 


351 


746 


Avgoustoglou et al. Ill] 




133770 




135196 








Savukov et al. [12] 




132738 




134231 








Dong et al. [13] 


134110 


134567 134940 


135969 


1402 


373 


830 


Froese Fischer and Tachiev [5] 134038 


134452 134807 


135887 


1435 


355 


769 


NIST [351 


134042 


134459 134819 


135889 


1430 


360 


777 



It should be pointed out, however, that the core excitations have been neglected in these 
two calculations. Moreover, in the work of Froese Fischer and Tachiev relativistic effects 
were included through the Breit-Pauli Hamiltonian but the orbit-orbit interaction, which is 
part of the Breit interaction, is ignored. In addition, we noticed that present calculated fine 
structure splittings are consistent with the experimental values [35|, and are better than 
other calculations as well. 
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B. The 31 Pf - El transitions 



In table |V] we report line strengths (S) and corresponding oscillator strengths (gf) in 
Babushkin and Coulomb gauges for 3 > l P° — 1 Sq transitions. These two gauges are related 
to the nonrelativistic length and velocity form of transition operators, respectively (37)]. 
The convergence of line strengths and oscillator strengths and the good consistency found 
between the two gauges further justify our computational models and suggest reliable atomic 
wave functions. 



Theoretical and experiments values published during the last two decades are also dis- 
played in table |V] For the 3 P° — 1 S'o transition, we see an excellent agreement with the 
semi-empirical calculations of Hibbert et al. Q and of Seaton [4]. The present gf value 
differs from MBPT values of Avgoustoglou et al. and of Savukov et al. 12] by 30% and 
18%, respectively. Such large discrepancies might be attributed to the Breit interaction that 
was completely or partly neglected in MBPT calculations of transition properties. Good 



13|. They 



agreement is found with the results obtained by Dong et al. 
approach 38j to account for non-orthogonal orbitals in transitions 



39 



adopted Lowdin's 



instead of the 



biorthogonal transformation technique used in this work. The difference between the results 

fl 

of Froese Fischer and Tachiev [5] and ours is about 12%. Using the similar Breit-Pauli 
approximation to Froese Fischer and Tachiev, Zatsarinny and Bartschat recently calculated 
the gf values by the B-spline method jl4], whose results approach our calculations. Com- 
pared with experimental measurements, our results perfectly agree with Zhong et al. and 
are in good consistency with Chan et al. 15| and Suzuki et al. 41] with respect to the 
experimental errors. It is worth noting that all these experiment measurements in good 
agreement with present calculations were obtained by the electron-energy-loss spectrometer 
(EELS) method. 



For the 1 P° — 1 Sq transition, the agreement between theories and experiments is better 
than for the spin-forbidden transition. But we find that the semi-empirical results of Hibbert 

4] and the B-spline values by Zatsarinny and Bartschat [l^] are larger 



et al. 



2| and of Seaton 



than other theoretical data. Present gf is also consistent with all experimental results listed 



in this table except for the value of Curtis et al. 



42]. 
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TABLE V. The convergence trends of line strengths S (in a.u.) and corresponding oscillator 
strengths gf for the 1,3 P° — 1 Sq El transitions of neutral neon. B: Babushkin gauge, C: Coulomb 
gauge. 



3 po 



So 



l p;- 1 s 



gf 



gf 



Model 


B 


C 


B 


C 


B 


C 


B 


C 


DF 


0.03172 


0.02584 


0.01360 


0.01108 


0.3428 


0.2796 


0.1484 


0.1211 


3SD 


0.03199 


0.03162 


0.01347 


0.01331 


0.3528 


0.3484 


0.1500 


0.1482 


4SD 


0.03250 


0.03315 


0.01325 


0.01351 


0.3474 


0.3539 


0.1431 


0.1458 


5SD 


0.03327 


0.03352 


0.01363 


0.01373 


0.3469 


0.3487 


0.1435 


0.1443 


6SD 


0.03383 


0.03378 


0.01387 0.01385 


0.3485 


0.3469 


0.1443 


0.1437 


7fSD 


0.03393 


0.03394 


0.01390 


0.01391 


0.3509 


0.3501 


0.1453 


0.1449 


8fSD 


0.03415 


0.03419 


0.01399 


0.01401 


0.3552 


0.3546 


0.1470 


0.1468 


9fSD 


0.03417 


0.03419 


0.01400 


0.01400 


0.3562 


0.3556 


0.1474 


0.1472 


CC 


0.03488 


0.03410 


0.01435 


0.01402 


0.3524 


0.3440 


0.1464 


0.1429 


MR 


0.03579 


0.03557 


0.01465 


0.01456 


0.3527 


0.3504 


0.1459 


0.1449 


Breit 


0.03032 


0.03007 


0.01241 


0.01231 


0.3583 


0.3556 


0.1482 


0.1471 



Hibbert et al. [2] 
Seaton [4] 
Avgoustoglou et al. [11] 
Savukov et al. [l£] 
Dong et al [13] 
Froese Fischer and Tachiev [5] 
Zatsarinny and Bartschat [14] 



Theories 

0.0123 

0.0126 

0.0163 0.0156 
0.0102 

0.03175 0.03309 0.01298 0.01353 
0.02680 0.01095 

0.0118 0.0116 



0.1607 

0.168 
0.161 0.147 
0.1459 

0.3492 0.3587 0.1442 0.1482 
0.3668 0.1514 

0.159 0.156 



Chan et al. [15] 
Ligtcnberg et al. [43J 
Suzuki et al. [41] 
Curtis et al. [42J 
Gibason et al. [44] 
Zhong et al. [451 



Experiments 

0.0118(6) 
0.01017(30) 

0.0106(14) 

0.0084(3) 
0.01095(32) 

0.0124(38) 



0.159(8) 
0.1369(35) 

0.137(18) 

0.165(11) 
0.1432(38) 

0.156(9) 



C. The 3 P 2 ° - 1 S M2 transition 

In table IVIl we display the ^P^ — l S$ M2 transition rates and corresponding line strengths 
as functions of the computational models as well as other theoretical and experimental values 
when available. It is found that our results are in good consistency with results of Beck jl9 ], 
Dong et al. [13j, Desclaux et al. [46] and Froese Fischer and Tachiev [5[. However, all 



theoretical predictions differ from the experimental value jl7j by amounts ranging from 14% 



to 40%. To explain such large discrepancies, further experiments are called for. 
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TABLE VI. Line strengths S (in a.u.) and rates A (in s _1 ) for the 3 P£ — 1 5'o M2 transition 
as a function of the active space. Numbers in square brackets stand for the power of 10, and in 
parentheses for error bars. 



Model 

1V1UUC1 


S 


A 




1/1 


3 766 


6 1QQ 


-21 
Z J 


*J -J 1 J 


3 730 


5 602 


-21 

Z J 


4SD 


3 Q1 6 


5 006 


-21 

Z J 


5SD 


4 031 


5 275 


-21 
Z J 


6SD 


4.159 


5.473 


-2] 


7fSD 


4.228 


5.548 


-2] 


8fSD 


4.332 


5.678 


-2] 


9fSD 


4.350 


5.697 


-2] 


CC 


4.284 


5.727 


-2] 


MR 


4.335 


5.657 


-2] 


Breit 


4.345 


5.672 


-2] 


mi 

1 heones 








Small- Warren and Chow Chiu [181 




4.10[- 


2] 


Indelicato et al. ' 




4.55[- 


2] 


Beck [12] 




5.05[- 


2] 


Dong et al. [131 




5.29[- 


2] 


Desclaux et al. * 




5.29[- 


2] 


Froese Fischer and Tachiev [5] 


4.525 


5.838 


-2] 



Experiments 



Zinner et al. [17] ^).06790(64) 
1" This values is cited in Ref. 47]; 
•f This values is cited in Ref. 17]. 

D. The 3 P ° - 3 Pf Ml and 3 P ° - 3 P 2 ° E2 transitions 

Line strengths and rates for 3 Pq — 3 P° Ml and 3 Pq — 3 P 2 ° E2 transitions are presented 
in table IVHI with the corresponding transition energies. For the Ml transition, we note 

12 



that the rate is much more sensitive to the transition energy than to the line strength that 
hardly changes with the computational models. As a result, higher-order electron correlation 
and the Breit interaction must be taken into account to achieve high accuracy for the Ml 
transition rate due to their considerably effects on fine structures as discussed in Sec. Ill Bl 
and Sec. Ill CI It is found from table IVHI that our final result is in good agreement with 
other theoretical calculations. 

For the E2 transition the rate is five orders of magnitude smaller than the Ml transition, 
and thus is negligible. However, we discovered that the transition probabilities in Babushkin 
and Coulomb gauges are not consistent with each other even with large configuration spaces. 
As can be seen from table IVTTj the inconsistency arises from the deviation of line strengths in 
Coulomb gauge from those in Babushkin gauge, although they converge with the expansion 
of configuration space. A strongly gauge-dependency of transition probabilities has also been 
found in the preceding investigation on the spin-forbidden 2s2p 3 P° — 2s 2 1 Sq transition 



of the Be-like C ion 



36 



48| . Chen et al. explained that this gauge dependency is caused 



by the neglect of the negative-energy state which significantly influence the velocity-gauge 



results 



49j. Therefore, we argue that the gauge dependency of the E2 transition rate in the 



case of Ne is brought about for the same reason. 



E. Hyperfine induced 3 Pq 2 ~ 1 Sq El transitions 

In the presence of hyperfine interactions, the electronic angular momentum J is coupled 
with the nuclear angular momentum / to form the total angular momentum F of the atomic 
system and only the latter is the good quantum number. As a result, new decay channels 
can be opened by hyperfine interactions, which affect lifetimes of metastable states substan- 
tially. These transitions, named as hyperfine induced transitions, have been investigated 
extensively during the last decade owing to their potential applications in many fields 20- 



25] . Neon possesses a stable isotope 21 Ne with nuclear spin / = 3/2, a magnetic dipole 
moment fii=— 0.661797 n.m. and with an electric quadrupole moment Q=0.103 barns in 



the nuclear ground state |50j]. Two El transitions from the metastable states 3 Pq 2 to the 
ground state l So can be induced by hyperfine interactions in 21 Ne isotope. In this section, 
we predict the decay rates of these two transitions. 



Methods calculating the HIT rate have been reviewed in Ref. 
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20] . Based on perturbation 



energies (in cm 1 ) of the 3 Pq 



3 i? Ml and 3 P ° 



TABLE VII. Line strengths S (in a.u.) and rates A (in s -1 ) together with corresponding transition 

3 Pg E2 transitions for neon. AE repre- 
sents transition energy. B: Babushkin gauge; C: Coulomb gauge. The number in square brackets 
represents the power of 10. 



Model 




Ml 






E2 






AE 


S A 


AE 




Sn 


An 


An 


DF 


389 


1 835 9 91 71- Si 


821 


S Q1 Ml 


1.26 


J- . \) >. ) 1 o 


5 271-81 

\J. — II O 


3SD 


387 


1 838 2 871 1-31 


806 


4.21[-1] 


1.61[-1] 


1.60 [-8 


6.14[-9] 


4SD 


449 


1 833 4 462 1-31 


854 


4.10[-1] 


6.61[-2] 


2.09[-8 


3.37[-9] 


5SD 


468 


1.829 5. 049 [-3] 


869 


4.12[-1] 


1.28[-1] 


2. 28 [-8 


7. 12 [-9] 


6SD 


474 


1.828 5.258[-3] 


875 


3.90[-l] 


1.65[-1] 


2. 24 [-8 


9.48 [-9] 


7fSD 


475 


1.828 5.293[-3] 


876 


3.72[-l] 


1.67[-1] 


2.15[-8 


9.61[-9] 


8fSD 


475 


1.830 5.279[-3] 


874 


3.20[-l] 


2.94[-2] 


1.82 [-8 


1.68 [-9] 


9fSD 


475 


1.830 5.296[-3] 


875 


3.09[-l] 


1.14[-2] 


1.77 [-8 


6.53[-10] 

L J 


cc 


479 


1.825 5.414[-3] 


880 


3.14[-1] 


3.83[-4] 


1.86 [-8 


2.27[-ll] 


MR 


408 


1.821 3.348[-3] 


845 


3.15[-1] 


1.47[-3] 


1.51[-8 


7. 08 [-11] 


Breit 


362 


1.849 2. 358 [-3] 


771 


3.14[-1] 


1.90 [-3] 


9.61 [-9 


5.81[-11] 


NIST 


359 




777 














Theory 












Small- Warren et al. [18] 




2.33[-3] 












Dong et al. [13] 




2.308[-3] 













Froese Fischer and Tachiev [5] 355 1.864 2. 240 [-3] 



theory, the HIT rate of Ne can be estimated by 

. 2.02613 x 10 18 a 

A = OHIT, (5) 

where A is the HIT transition wavelength in A, Shit the corresponding line strength that 
is expressed as 

5Wjr=IM s ^l|0 (1) || 1 5'o> + /i a < 1 ^l|0 (1) ir5 , o>| a . (6) 

For the latter equation, we only take into account the effect of the adjacent 3 P° and 1 P° 
perturbative states. The two reduced matrix elements appearing in Eq. (jSJ) are the square 
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TABLE VIII. Hyperfine induced 3 Pq — 1 Sq El transition rates A (in s _1 ) for 21 Ne together with 
off-diagonal hyperfine interaction matrix elements W in (a.u.) and hyperfine mixing coefficients as 
functions of computational models. The number in square brackets represents the power of 10. 



( 3 Pf , 3 P °) 3 P °) 



Model 




W\ 




fix 






VV2 








A 

A 


i \ i ■ 
Ub 


— 1 


/I O A 1 

4z41 


"'J 


—8.1)32 


-5J 


— 1 


zzyu 


-7J 


2.668 


-5J 


1.716 


SbD 


— 1 


3692 


"'J 


— 7.767 


rl 


— 1 


262b 


-7J 


2.730 


-5J 


i a A A 

1.644 


4bD 


— 1 


1UU1 


"'J 


—5.383 


Kl 

-5J 


— 1 


onon 
oUol) 


-7J 


O.U28 


rl 

-oj 


l.ZOZ 


^F) 

oou 


1 

JL 


-LOUD 


71 


— 5 497 


51 


_1 


9&QQ 


71 


O.UOi 


51 


1 ^1 fi 
1 .oxu 


6SD 


-1 


1210 


-7] 


-5.188 


-5] 


-1 


2951 


-7] 


3.140 


-5] 


1.323 


7fSD 


-1 


1343 


-7] 


-5.238 


-5] 


-1 


2918 


-7] 


3.124 


-5] 


1.328 


8fSD 


-1 


1288 


-7] 


-5.219 


-5] 


-1 


2870 


-7] 


3.079 


-5] 


1.312 


9fSD 


-1 


1297 


-7] 


-5.218 


-5] 


-1 


2878 


-7] 


3.083 


-5] 


1.316 


CC 


-1 


3475 


-7] 


-6.173 


-5] 


-1 


2337 


-7] 


3.000 


-5] 


1.458 


MR 


-1 


3162 


-7] 


-7.072 


-5] 


-1 


2411 


-7] 


2.776 


-5] 


1.488 


Breit 


-1 


3438 


-7] 


-8.156 


-5] 


-1 


2072 


-7] 


2.614 


-5] 


1.484 



roots of line strength S presented in table [V] h\ and h 2 in Eq. (jBJ) stand for the hyper- 
fine mixing coefficient that can be estimated from the ratio of the off-diagonal hyperfine 
interaction matrix element and the energy difference between the interactive states. 

Using the computational model described in Sec. Ill B[ we calculate the hyperfine in- 
duced 3 Pq 2 — 1 Sq El transition rates and present the results in table IVIIII and table IIXI 
Additionally, the off-diagonal hyperfine interaction matrix elements {W) and the hyperfine 
mixing coefficients are displayed as well. It is found from table IVIIII that the off-diagonal 
hyperfine interaction matrix elements are well converged with the expansion of the config- 
uration space. While relatively large changes in the hyperfine mixing coefficients between 
CC, MR and Breit models are mainly attributed to the energy separations involved that 
are sensitive to the higher-order correlation effects as discussed in Sec. Ill Bl As can be seen, 
the final hyperfine induced transition rate is three orders of magnitude larger than the Ml 
transition presented in Sec. IIII Dl and thus reduces the lifetime of the states by a factor of 
630. Therefore, for 21 Ne the HIT is a dominant decay channel from the 3 Pq state. 
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TABLE IX. F-dependent hyperfine induced — 1 Sq transition rates A (in s _1 ) together with 
associated hyperfine mixing coefficients hi and I12 for 21 Ne by using the "Breit" model. The number 
in square brackets represents the power of 10. 

F hi h 2 A 

1/2 -4.946[-7] 5.089[-7] 2.500[-4] 

3/2 5.779[-6] 1.351[-6] 6.395 [-5] 

5/2 1.935[-5] 2.390[-6] 6.153[-3] 

For the other hyperfine induced transition from the 3 P 2 ° state to the ground state, 
the mechanism is a little more complex since the excited level possesses several hyper- 
fine sublevels with F = 1/2,3/2,5/2,7/2 for the 21 Ne isotope. Out of them only the 
F = 1/2,3/2,5/2 states can decay to the ground state. In table [IX] we present the transi- 
tion rates and corresponding hyperfine mixing coefficients for these hyperfine states using 
the Breit model. As can be seen from this table, the HIT rates are somewhat smaller than 
the M2 transition probability discussed in Sec. IIII CI but still significantly affect the level 
lifetime. 



F. Level lifetimes in 2s 2 2p 5 3s configuration 

Using the data presented in tables |V] - [IX] we obtain the lifetimes of states in 2p 5 3s 
configuration for 20 ' 21 Ne isotopes by 

where the summation is made over the main decay channels. For the 3 P| state of 21 Ne 
isotope, the weighted average lifetime (r = ^^|q^y ) is calculated. The results are reported 
in table |2D It can be seen that the lifetimes of those two metastable states are apparently 
different owing to the impact of hyperfine interactions, especially for the 3 Pq state. We 
should emphasize that the interference effect between the main decay channels is neglected 
in Eq. (J7J), which brings about observable variation in lifetimes if transition probabilities 
have similar orders of magnitude. As discussed in Sec. IIII El the hyperfine induced transition 
rate of the 3 P 2 ° state for 21 Ne has the same order of magnitude as the M2 transition, and 
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TABLE X. Lifetimes (in s) of levels in 2p 5 3s configuration for 20,21 Ne isotopes. The relevant nuclear 

□ 

parameters are taken from Ref. [50J. The number in square brackets represents the power of 10. 



isotope 


3 no 3 no 3 po 
r 2 r \ M) 


Ipo 


20 Ne 


17.63 1.995[-8] 424.1 


1.638 [-9] 


21 Ne 


17.10 1.995[-8] 0.6728 


1.638 [-9] 



strong interference may occur. This also influences the radiative emission distribution, which 



is useful for anisotropy plasma diagnosis 



5l|. Further studies are ongoing. 



G. Estimation of uncertainties 

For light atoms such as neon the main uncertainties in calculations of physical quantities 
arise from electron correlation effects. In this work, large-scale configuration spaces are used 
to account for these correlation effects in the case of neutral neon, even partly including 
higher-order correlation among 2s, 2p valence electrons. The residual higher-order valence 
correlations and the higher-order correlations between Is core electrons and between core 
and valence electrons, which are not taken into account, contribute to the uncertainties. 
By monitoring the convergence of physical quantities under investigation as the active set is 
enlarged as well as monitoring the changes as the correlation models are defined by including 
higher-order correlation effects, we estimate that the errors in present calculations is about 
2%. This observation is further strengthened by the excellent agreement between El rates 
in the length and velocity gauges. The hyperfine induced 3 P| — 1 S'o El transition rate is 
an exception. This transition is sensitive to higher-order correlation effects not included or 
saturated in our calculations. Moreover, the counteraction between off-diagonal magnetic 
dipole and electric quadrupole interactions contributes to the uncertainties in this rate. 
Approximately, these bring about 10%-20% error for this transition rate. Other physical 
effects neglected in this work such as frequency-dependent Breit interactions and quantum 



electrodynamical (Q 
Avgoustoglou et al. 



3D) corrections are indeed fractional for neutral neon, as discussed by 

10]. 
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IV. CONCLUSION 



In this work we investigate the transition properties of the main one-photon decay chan- 
nels for the 2p 5 3s configuration of Ne isotopes using the MCDHF method. The electron 
correlation effects are taken into account systematically with the active space approach. 
Detailed comparisons are made with measurements and with other calculations. The effects 
of Breit interaction on fine structures and transition properties are discussed. It is found 
that the Breit interaction changes the line strength of the 3 P° — 1 Sq transition by around 
17%. Present calculations do not resolve the discrepancies in the 3 P 2 ° — 1 S M2 transi- 
tion rates between theories and experiments. Further measurement is therefore called for. 
The hyperfine induced 3 P ° 2 — 1 S El transition rates for the 21 Ne isotope are calculated 
as well. We discovered that the hyperfine interactions drastically affect the lifetime of the 
metastable states, especially for the 3 Pq state. The lifetime of states in 2p 5 3s configuration 
are predicted for both 20 Ne and 21 Ne isotopes. 

ACKNOWLEDGEMENT 

This work was supported by the Communaute frangaise of Belgium (Action de Recherche 
Concertee), the Beglian National Fund for Scientific Research (FRFC/IISN Convention) and 
by the IUAP Belgian State Science Policy (Brix network P7/12). PJ and GG acknowledge 
support from the Visby program of the Swedish Institute. 



18 



[1] D. L. Matthews, P. L. Hagelstein, M. D. Rosen, M. J. Eckart, N. M. Ceglio, A. U. Hazi, 
H. Medecki, B. J. Macgowan, J. E. Trebes, B. L. Whitten, E. M. Campbell, C. W. Hatcher, 
A. M. Hawryluk, R. L. Kauffman, L. D. Pleasance, G. Rambach, J. H. Scofield, G. Stone, 
and T. A. Weaver, Phys. Rev. Lett. 54, 110 (1985). 

[2] A. Hibbert, M. Le Dourneuf, and M. Mohan, Atomic Data and Nuclear Data Tables 53, 23 
(1993). 



[3] 
[4] 
[5 
[6 
[7 
[8 
[9 

[io; 
[ii 

[12 
[13 

[14 

[15 
[16 
[17 

[18 
[19 
[20 
[21 



A. Hibbert and M. P. Scott, J. Phys. B: At. Mol. Opt. Phys. 27, 1315 (1994). 
M. J. Seaton, J. Phys. B: At. Mol. Opt. Phys. 31, 5315 (1998). 



C. Froese Fischer and G. Tachiev, Atomic Data and Nuclear Data Tables 87, 1 (2004) . 



T. Morel and K. Butler, A&A 487, 307 (2008) 



K. Butler, A&G 49, 6.23 (2008). 

W. Fielder, Jr., D. L. Lin, and D. Ton-That, Phys. Rev. A 19, 741 (1979). 

W. Fielder, Jr., D. L. Lin, and D. Ton-That, Phys. Rev. A 20, 1738 (1979). 

E. Avgoustoglou, W. R. Johnson, Z. W. Liu, and J. Sapirstein, 



Physical Review A 51, 1196 (1995) 



E. N. Avgoustoglou and D. R. Beck, Physical Review A 57, 4286 (1998) 



I. M. Savukov, W. R. Johnson, and H. G. Berry, Physical Review A 6 6, 052501 (2002) 

C. Z. Dong, L. Y. Xie, X. X. Zhou, X. W. Ma, and S. Fritzsche, Hyperfine Interactions 

146/147, 161 (2003). 



O. Zatsarinny and K. Bartschat, Phys. Scr. T134, 014020 (2009) 



W. F. Chan, G. Cooper, X. Guo, and C. E. Brion, Phys. Rev. A 45, 1420 (1992). 



R. S. Van Dyck, Jr., C. E. Johnson, and H. A. Shugart, Physical Review A 5, 991 (1972) 
M. Zinner, P. Spoden, T. Kraemer, G. Birkl, and W. Ertmer, 

|Physical Review A 67, 010501 (R) (2003) 



N. E. Small- Warren and L. Y. Chow Chiu, Phys. Rev. A 11, 1777 (1975). 



D. R. Beck, Phys. Rev. A 66, 034502 (2002) 



W. R. Johnson, Can. J. Phys. 89, 429 (2011) 



P. Indelicato, F. Parente, and R. Marrus, Phys. Rev. A 40, 3505 (1989). 

19 



[22] M. Andersson, K. Yao, R. Hutton, Y. Zou, C. Y. Chen, and T. Brage, 

|Phys. Rev. A 77, 042509 (2008) [ 
[23] J. G. Li, P. Jonsson, G. Gaigalas, and C. Z. Dong, 



The European Physical Journal D 51, 313 (2009) 



[24] J. G. Li and C. Z. Dong, Plasma Science and Technology 12, 364 (2010). 
[25] M. H. Chen and K. T. Cheng, |Can. J. Phys. 482, 473 (2011)[ 



[26] P. Jonsson, X. He, C. Froese Fischer, and I. P. Grant, 

|Computer Physics Communications 177, 59 7 (2007) , 
[27] I. P. Grant, Relativistic Quantum Theory of Atoms and Molecules (New York: Springer, 2007). 
[28] J. Olsen, M. R. Godefroid, P. Jonsson, P. A. Malmqvist, and C. Froese Fischer, Phys. Rev. 

E 52, 4499 (1995). 

[29] S. Verdebout, P. Jonsson, G. Gaigalas, M. Godefroid, and C. Froese Fischer, 
| J. Phys. B: At. Mol. Opt. Phys. 43, 074017 (20ToJ 



[30] P. Jonsson, T. Carette, M. Nemouchi, and M. Godefroid, 

|J. Phys. B: At. Mol. Opt. Phys. 43, 115006 (2010)1 



[31] T. Carette and M. R. Godefroid, J. Phys. B: At. Mol. Opt. Phys. 44, 105001 (2011) 



[32] T. Carette and M. R. Godefroid, Phys. Rev. A 83, 062505 (2011) 



[33] C. Froese Fischer, T. Brage, and P. Jonsson, Computational Atomic Structure An MCHF 

Approach (London: Institute of Physics Publishing, 1997) p. 290. 
[34] T. Lee, N. C. Dutta, and T. P. Das, Phys. Rev. A 4, 1410 (1971). 

[35] A. Kramida, Y. Ralchenko, J. Reader, and N. A. Team, 



"NIST Atomic Spectra Database (version 5.0) [Online]," (2012). 



[36] A. Ynnerman and C. Froese Fischer, Phys. Rev. A 51, 2020 (1995). 
[37] I. P. Grant, J. Phys. B: Atom. Molec. Phys. 7, 1458 (1974). 
[38] P.-O. L6wdin, Phys. Rev. 97, 1490 (1955). 

[39] S. Fritzsche, C. Froese Fischer, and C. Z. Dong, Computer Physics Communications 124, 
340 (2000). 



[40] S. Fritzsche, Computer Physics Communications 183, 1525 (2012) 



[41] T. Y. Suzuki, H. Suzuki, S. Ohtani, B. S. Min, T. Takayanigi, and K. Wakiya, Phys. Rev. A 
49, 4578 (1994). 



20 



[42] L. J. Curtis, S. T. Maniak, R. W. Ghrist, R. W. Irving, D. G. Ellis, M. Henderson, M. H. 
Kacher, E. Trabert, J. Granzow, P. Bengtsson, and L. Engstrom, Phys. Rev. A 51, 4575 
(1995). 

[43] R. C. G. Ligtenberg, P. J. M. van der Burgt, S. P. Renwick, W. B. Westerveld, and J. S. 

Risley, Phys. Rev. A 49, 2363 (1994). 
[44] N. D. Gibson and J. S. Risley, Phys. Rev. A 52, 4451 (1995). 

[45] Z. P. Zhong, S. L. Wu, R. F. Feng, B. X. Yang, Q. Ji, K. Z. Xu, Y. Zou, and J.-M. Li, Physical 

Review A 55, 3388 (1997). 
[46] This value is cited in Ref. 

[47] M. Walhout, A. Witte, and S. L. Rolston, Phys. Rev. Lett. 72, 2843 (1994). 

[48] P. Jonsson and C. Froese Fischer, Physical Review A 57, 4967 (1998). 

[49] M. H. Chen, K. T. Cheng, and W. R. Johnson, [Physical Review A 64, 042507 (2001) 

[50] N. J. Stone, |Atomic Data and Nuclear Data Tables 90, 75 (2005) [ 

[51] K. Yao, M. Andersson, T. Brage, R. Hutton, P. Jonsson, and Y. M. Zou, 
Phys. Rev. Lett. 97, 183001 (200"6)| 



21 



